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ABSTRACT 

Context. Miras are long-period variables thought to be in the asymptotic giant branch (AGB) phase of evolution. In about one percent 
of known Miras, the pulsation period is changing. It has been speculated that this changing period is the consequence of a recent 
thermal pulse in these stars. 

Aims. We aim to clarify the evolutionary state of these stars, and to determine in particular whether or not they are in the thermally- 
pulsing (TP-)AGB phase. 

Methods. One important piece of information that has been neglected so far when determining the evolutionary state is the presence 
of the radio-active s-process element technetium (Tc). We obtained high-resolution, high signal-to-noise-ratio optical spectra of a 
dozen prominent Mira variables with changing pulsation period to search for this indicator of TPs and dredge-up. We also use the 
spectra to measure lithium (Li) abundances. Furthermore, we establish the evolutionary states of our sample stars by means of their 
present-day periods and luminosities. 

Results. Among the twelve sample stars observed in this programme, five were found to show absorption lines of Tc. BH Cru is 
found to be a carbon-star, its period increase in the past decades possibly having stopped by now. We report a possible switch in the 
pulsation mode of T UMi from Mira-like to semi-regular variability in the past two years. R Nor, on the other hand, is probably a 
fairly massive AGB star, which could be true for all meandering Miras. Finally, we assign RU Vul to the metal-poor thick disk with 
properties very similar to the short-period, metal-poor Miras. 

Conclusions. We conclude that there is no clear correlation between period change class and Tc presence. The stars that are most 
likely to have experienced a recent TP are BH Cru and R Hya, although their rates of period change are quite different. 

Key words. Stars: AGB and post-AGB - Stars: oscillations ~ Stars: evolution - Stars: abundances 



1. Introduction 

Mira-type variables are red giant stars that exhibit long-period 
variability with large amplitude (AV > 2™5). The pulsation pe- 
riod of most Mira variables is stable over long time spans (i.e. of 
between decades and centuries), but there are exceptions. |Zijlstra| 
|& Bedding (2002l define three classes of secular period change 
in Mira variables: 

Continuous change: A continuous decrease or increase in the 
period over a century or more, with no indication for epochs 
of a stable period. The change in period is in all cases large, 
15% or more over 100 years. In the following, we refer to 
this class as CCh. 

Sudden change: After a long phase of stability, the period sud- 
denly begins to either decrease or increase at a large rate. 



* Based on observations made with the Mercator Telescope, operated 
on the island of La Palma by the Flemish Community, at the Spanish 
Observatorio del Roque de los Muchachos of the Instituto de Astrofisica 
de Canarias, and with the Swiss 1.2 m Euler telescope at La Silla, Chile. 
The spectra are only available in electronic form at the CDS via anony- 
mous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via |http://cdsweb.i>l 
|strasbg.fr/cgi-bin/qcat?J/A-HA/| 



The total change in period is similar to that of the CCh class, 
but the rate is ten times faster because the same change is 
achieved within a decade rather than within a century (SCh 
in the following, not to be confused with the SC spectral 
type!). 

Meandering change: A number of long-period Miras show ev- 
idence of meandeiing or fluctuating periods. The periods 
change by 10% over several decades, followed by a return to 
the previous period. The rate of change is comparable to that 
of the CCh Miras, but the total change is somewhat smaller 
than that of the other two classes. Most Miras in this class 
have periods in excess of 400 d (MCh in the following). 



Templeton et al. ( 2005 1 study the period evolution of 547 



Mira variables, of which 57 are found to have a period change at 
the 2cr confidence level, and eight stars at the level of 6cr. Their 
catalogue is the current reference for observed secular evolution 
in Miras in the solar neighbourhood. 

From the point of view of stellar evolution, Mira vaiiables 
are thought to be in the asymptotic giant branch (AGB) phase. 
The hydrogen-burning shell is the main energy source dur- 
ing most of the time on the AGB, which is quasi-periodically 
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interrupted by powerful ignitions of the helium-burning shell 
( Schwarzschild & Harm 1965]). These ignitions are called 



helium-shell flashes or thermal pulses (TPs), hence this phase is 
also referred to as the TP-AGB. A TP has a large impact on the 
structure of the whole star: the luminosity, temperature, radius, 
and thus pulsation period of the outer envelope, are predicted to 
vary strongly within centuries. For a star with 2 M© in the middle 
of its TP-AGB phase, the TP itself (convective zone in the inter- 
shell region) lasts about 300 - 400 years. The evolution during 
the TP is governed by the thermal timescale of the envelope, and 
the evolution is faster in both later TPs and more massive stars. 
Several hundred years after the onset of the TP, surface abun- 
dance changes can occur as the third dredge-up (3DUP) begins 
to operate. This mixing event may take place after a TP when 
the convective envelope deepens to reach layers that have pre- 
viously undergone nuclear processing ( |Busso et al.||T999l ). The 
onset of 3DUP is predicted to coincide with a temporary lumi- 
nosity maximum. Once the convective envelope reaches layers 
with processed material, it should become visible on the stellar 
surface very quickly because the turn-over time of the convective 
envelope is of the order of only a few years. 

The first suggestion that the observed period change in a few 
Mira stars could be related to a recent TP dates back to IWoodI 
( |1975| l. [Wood & Zarro| ( |1981| l present evolutionary models of 
the TP-AGB and attempt to derive constraints on the core mass 
and luminosity of three Miras (R Hya, R Aql, and W Dra) with a 
well-known, long-term period change. The luminosity evolution 
inferred from the observed period change of the stars is repro- 
duced well by the predicted luminosity evolution of the models. 
This gives some support to the suggestion that a TP is the cause 
of the period change. The fraction of Miras with large period 
changes (~1.6%; Templeton et al. 20051 agrees well with the 
expected ratio of pulse duration to inter-pulse time of roughly 1 - 
2%. Hence, from a statistical point of view, we expect that at any 
given time at least 1 % of the Miras is undergoing a TP. Whether 
the ones that are observed to show a large period change are 
also the ones that currently undergo a TP is a different question, 
which we aim to address here. 

There are, however, also alternative explanations of the ob- 



served period changes, which do not require a recent TP. Ya'ari 
|& TuchmanI ([T9%| and [Lebzelter & Wood] ( |2005] l identify a 
feedback mechanism between the pulsation and the stellar en- 
tropy structure, which can lead to mode switching. On the other 
hand, [Z^lstra et al] ( |2004[ ) propose that a possibly chaotic feed- 
back between molecular opacities, pulsation amplitude, and pe- 
riod can cause an unstable period. This is supposed to be most 
important for stars with a C/0 ratio very close to unity (SC and 
CS spectral types), as small changes in the temperature in the at- 
mosphere can cause large changes in the molecular abundances, 
and hence the opacities. This alternative explanation might be 
applicable to the period change in the MCh group, as their peri- 
ods change back and forth on a much shorter timescale than can 
be explained by a TP. 

A piece of information that has not yet been taken into ac- 
count in this discussion is the presence of the radio-active ele- 
ment technetium (Tc) in the atmospheres of these stars. Tc is 
a product of the slow neutron-capture process (s-process) that 
takes place in deep layers of AGB stars. The s-process prod- 
ucts are subsequently mixed to the surface by a 3DUP event. 
The longest-lived Tc isotope produced by the s-process is "^^Tc 
with a half-life time of 210000 years, which is short compared 
to stellar evolution timescales before the AGB. If observed on 
the surface of a star, Tc is an indisputable sign of an ongoing 
s-process and 3DUP in that star, hence indicates that the star is 



indeed in the TP-AGB phase of evolution. Furthermore, 3DUP 
mixes up carbon, the product of He-burning, which can change 
the star from an oxygen-rich M-type to a carbon-rich C-type star. 
In previous surveys of Tc in red giant stars (e.g. [Little et| 



al. 1987) 1, only a few Miras with changing pulsation periods 



were investigated, and for those that had been observed the pres- 
ence of Tc was often unclear. The presence of Tc would be sup- 
port for the hypothesis of a recent TP that might be the cause 
of the observed period change. We note, however, that the ab- 
sence of Tc does not exclude that TPs are going on in a given 
star Nevertheless, the absence of Tc places a constraint on the 
strength of the TP, because only the strongest TPs on the upper 
AGB will be followed by a 3DUP event. 

In addition to Tc, the light element lithium (Li) is an impor- 
tant diagnostic tool for stellar evolution studies, in particular for 
internal mixing processes. It is very fragile because it is burnt to 
helium by proton captures at temperatures exceeding 3 x 10^ K. 
Standard stellar evolution predicts that Li should diminish dur- 
ing the evolution to the giant phase, and indeed red giants are 
generally observed to have a low Li abundance compared to 
main sequence stars. However, standard stellar evolution also 
predicts that in high mass AGB stars (M > 4-Mq), Li may be 
produced in large amounts. In these stars, the bottom of the con- 
vective envelope extends to the H-burning shell (also called hot 
bottom burning, or HBB; Sackmann & Boothroyd I992[l, which 
activates the Cameron-Fowler mechanism ( [Cameron & Fowler| 
1971 ) producing Li. Moreover, these stars do not become C-rich 
by means of the dredge-up of C, because the CN cycle also con- 
verts '^C to '"^N. Li surface abundances of up to log e(Li) - 4.5 
(where loge(Li) - log A^(Li)/A^(H) + 12) are predicted (e.g. 
Sackmann & Boothroyd 1992 Karakas 2003 1. Luminous, Li- 
rich stars have been observed in both the Magellanic Clouds 
(Smith et al. 1995) and the Galaxy ( Garcfa-Hemandez et al. 



2007) . Hence, the observation of a luminous O- and Li-rich star 
infers a fairly high mass for that star. 

Lower-mass AGB stars can also be somewhat enriched in Li. 
The layer of radiative energy transport, which usually prevents 
the operation of the Cameron-Fowler mechanism in these stars, 
is naturally overcome during 3DUP events and can lead to Li 
enhancement at the surface up to log e(Li) < 1.8 ( [Karakas et al. 



2010 1. In addition to this, non-convective extra-mixing processes 
(e.g. thermohaline mixing, magnetic buoyancy) may provide 
transport of matter across the radiative barrier. Observational ev- 



idence along these lines is presented by Vanture et al 
Utten flialer et a l.[( [2007bl l, and [Uttenflialer & Lebzelter 



(2007: 



(2010 1 



However, measurements of the Li abundances of O-rich, M-type 
AGB stars remain scarce. 

For the present work, we selected several prominent repre- 
sentatives of the three period change classes described above 
and observed their optical spectra at high resolution. With these 
data, we can safely establish the evolutionary state of the Tc- 
rich stars on the TP-AGB. A central aim of the paper is to check 
whether the presence of Tc is correlated with the different period 
change classes. The spectra also provided us with the opportu- 
nity to measure the abundance of Li. In addition, we determined 
the present-day pulsation periods of the stars from AAVSO ob- 
servations, to determine whether the period change has contin- 
ued in recent years. To construct an HR diagram of the stars, we 
calculated their luminosities with a period-magnitude relation. 

The structure of the paper is as follows. In Sect.|2j the target 
selection and observations are described; Sect. [3] is devoted to 
the analysis of the spectra and the results thereof; the results are 
discussed in Sect.|4] and finally Sect.|5]presents our conclusions. 
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Table 1. Observation log. 



Name 


Instr. 


Date 


Exp. time 


SNR 






(dd/mm/yyyy) 


(s) 


(at430nm) 


BH Cru 


C 


03/08/2009 


4x 1800 


12 


RU Vul 


H 


04/07/2009 


300 + 2 X 600 


30 


TUMi 


H 


04/07/2009 


4 X 1800 


20 


R Aql 


H 


04/07/2009 


900 + 2 X 300 


50 


RHya 


H 


13/02/2010 


2x400 


60 


WDra 


H 


15/08/2009 


3 X 1800 


~0 


RNor 


C 


17/08/2010 


2 X 1800 


20 


SHer 


H 


22/10/2009 


630 + 2 X 1200 


30 


SOri 


H 


13/02/2010 


3 X 1800 


15 


TCep 


H 


25/10/2009 


2 X 1200 


140 


WHya 


H 


13/02/2010 


3x 150 


50 


THer 


H 


04/07/2009 


2x 800 


70 


Note: The 


meaning 


of the columns 


is as follows: (1) Name listed in 



Simbad; (2) Instrument used (H: Hermes, C: Coralie); (3) Mean date 
of the observations; (4) Exposure times of the spectra in seconds; (5) 
Approximate SNR at 430 nm. 



2. Target selection and observations 

The targets of our observations were selected from the stars dis- 
cussed in [Wood & Zarol l pMT) and |Zijlstra & Bedding] ( [2002) l. 
They are prominent representatives of the period change classes 



introduced by |Zijlstra & Bedding (2002 ), and most are also dis- 
cussed in detail in |Templeton et al., ( ,2005 ) . In total, twelve stars 
were selected for the observations. 

The targets that are most readily observable from the north- 
ern hemisphere were observed with the Hermes spectrograph 
( [Raskin et al.| |2010| l, which is mounted on the 1.2 m Mercator 
telescope at the Roque de los Muchachos observatory on the is- 
land of La Palma, Spain. The two targets best observable from 
the southern hemisphere were observed with the Coralie spec- 
trograph, which is mounted on the 1.2 m Leonhard Euler tele- 
scope at La Silla observatory, Chile. Both Hermes and Coralie 
are fibre-fed, optical high-resolution spectrographs, mounted on 
twin telescopes. Hermes observes the range 377 - 900 nm at a 
spectral resolution of R - A/AA - 85 000, whereas Coralie ap- 
proximately covers the range 385 - 690 nm, with a somewhat 
lower spectral resolution of 7? = 50 000. In terms of transmis- 
sion, Hermes is by far the more efficient of the two instruments. 

The sample stars and the log of the observations is presented 
in Table [T] We aimed to observe the stars close to their visual 
maximum brightness across their pulsation cycle, to maximise 
the signal-to-noise ratio (SNR) of the spectra. 

The observed 2D spectra were reduced with the respective 
instrument pipeline software. As a second measure to increase 
the SNR in the blue spectral range, we combined the individual 
ID spectra into one. This was essential for some of the stars, as 
they have a very low flux in the blue range. Despite these efforts 
to increase the SNR, most stars have a low to moderate SNR at 
430 nm. W Dra is on average the faintest star in the sample and 
could not be observed close to its visual brightness maximum. 
As a result, the SNR of its spectrum in the B-band is essentially 
zero. Nevertheless, the signal in the red part of the spectrum was 
high enough to search for other signatures of s-process enrich- 
ment (see below). The SNR at 671 nm, where the Li resonance 
doublet lies, is about 25 in W Dra's combined spectrum, which 
is just good enough to determine the Li abundance. For the other 
stars, the SNR in this spectral range is very high (> 200), so we 
used one of the individual spectra instead of a combined spec- 
trum to measure the Li abundance (see Sect. |3.1.4] i. 



In Table [T] the stars are grouped into the different period 
change classes described in the introduction, and separated by 
horizontal lines: The first group is that of the sudden change stars 
(SCh), the second one that of the continuous change stars (CCh), 
and the third one the meandering Miras (MCh). The star T Her 
in the last row is an exception here. It is a Mira variable with 
a constant period, hence it does not belong to any of the period 
change classes. However, it is discussed among other stars in 
[Wood & Zarro| ( |T981|l, and it was hitherto uncertain whether this 
star contains Tc. Little et al. ( 1987 1 classify it as "Tc doubtful", 
and we wish to clarify this situation. Furthermore, RU Vul and 
W Hya are actually classified as SRVs in the Simbad database. 
This classification is controversial, as we will further elaborate 
on in Sect.lO 



3. Analysis and Results 

3.1. Analysis of the spectra 

3.1 .1 . Oxygen-rich versus carbon-rich 

As a first step, we performed a global inspection of the spec- 
tra to decide whether they are dominated by either features of 
oxygen-bearing molecules (e.g. TiO, VO) or those of carbon- 
bearing molecules (e.g. C2, CN). Stars whose spectra are domi- 
nated by TiO and VO bands are oxygen-rich (i.e. they have a C/O 
ratio < 1 in their atmosphere), whereas stars whose spectra are 
dominated by features of C2 and CN are carbon-rich (C/O > 1). 
It turns out that all but one star in our sample displays TiO bands, 
this one star being BH Cru. This star has prominent C2 and (less 
prominent) CN features in its spectrum, whereas ZrO bands are 
absent. Additional very prominent features in BH Cru include 
the Na D line, indicating a very low temperature and/or C/O ra- 
tio, and the Ball 4554 and Sri 4607 lines, indicating an s-process 
enrichment. Figure [T] shows the Coralie spectrum of BH Cru, 
smoothed to the resolution of classification spectra to make the 
main features visible. The C2 bands clearly dominate the spec- 
trum. 

Among the other stars, the one with clearly the weakest TiO 
features is RU Vul, implying that it is the hottest star in our 
sample (see also Sect. |3.1.4| . All other O-rich stars display pro- 
nounced bands of TiO. The range of spectral types found in the 
compilation of Skiff| (2010l is listed in column 3 of Table |2] All 
stars exhibit hydrogen emission lines from time to time. 

3.1.2. Technetium content 

Before searching for Tc, it is necessary to measure the radial 
velocity of the stars. To this end, we applied a cross-correlation 
technique, using a synthetic model spectrum as a template. The 
spectral range of 430 - 435 nm was used to complete this, ex- 
cept for W Dra, for which only the red part of the spectrum had 
enough signal and we then used the range 701 - 710 nm. The 
found heliocentric radial velocities are summarised in column 4 
of Table |2] 

With the radial velocity information, we determined whether 
or not Tc is present in the sample stars using the flux-ratio 
method introduced by [Uttenthaler et a l. (2007a), and by visual 
inspection of the spectra around the strongest lines of neutral 
Tc (A = 398.497, 403.163, 404.911, 409.567, 423.819, 426.227, 
and 429.706 nm, in air). In this method, the mean flux contained 
within a range of a few wavelength pixels centred on a quasi- 
continuum point is divided by the mean flux in a small region 
around the Tc line. If Tc line absorption is present, the result- 
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700 



Wavelength (nm) 



Fig. 1. Coralie spectrum of BH Cm, smoothed to a resolution of 
~0.3 nm per wavelength point. The main features of the spec- 
trum are identified. The inset shows a zoom on the spectrum 
blue-ward of 490 nm, where the flux is very low. 




1.0 1.5 

cont. / line Tc 426. 2nm 

Fig. 2. Continuum-to-line flux ratios of the Tc lines at 426.2 and 
429.7 nm of our sample stars. The data points of the Tc-rich 
stars, i.e. those stars with an enhanced ratio in both hnes, are 
marked with their names. 



ing ratio will be considerable larger than one, but close to one 
otherwise. The measured flux ratios for the two Tc lines at the 
longest wavelengths, hence highest SNR, are shown in Fig. [2] 
For the Tc 426.2 nm line, the range 426.120 - 426. 140 nm was 
used for the continuum point and 426.220 - 426.240 nm for the 
line, while for the Tc 429.7 nm line the ranges were 429.625 - 
429.655 and 429.690 - 429.718 nm, respectively. The error bars 
of the data points were determined with a Monte Carlo simula- 
tion. One-hundred trial spectra were created by adding Gaussian 
noise with the magnitude of the inverse SNR to the observed, 
combined spectrum. The flux ratios were measured from those 
trial spectra, and the standard deviation of the distribution of the 
ratios was adopted as the error bar. 

Five stars were found to have a higher continuum-to-line flux 
ratio, which was confirmed by visual inspection of the spectra 
around the Tc lines. These stars are BH Cm, R Hya, S Her, S Ori, 
and T Cep. The surprise in this list is R Hya, which was previ- 
ously classified by Uttenthaler & Lebzelter (2010 1 as Tc-poor 
on the basis of a UVES spectrum of that star. After careful in- 
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423.70 423.80 423.90 
Wavelength (nm) 



424.00 



Fig. 3. Spectra of three sample stars around the Tc line at 
423.8 19 nm. From top to bottom: T UMi, R Hya, and S Her. 
The spectra of T UMi and R Hya have been shifted upwards by 
1.2 and 0.6 to make the plot clearer The wavelength of the Tc 
line is marked with a dotted vertical line. 



spection of all Tc lines in the observed range and a comparison 
between the UVES and the Hermes spectrum of R Hya, it be- 
came clear that this star does indeed have Tc in its atmosphere 
in contrast to the findings of Uttenthaler & Lebzelter| ( |2010 



The UVES observations on 6 July 2000 were performed near 
maximum light of R Hya, while the Hermes observations were 
done about half-way between a maximum and a minimum. A 
cross-correlation of the two spectra shows that there are veloc- 
ity components in the atmosphere that differ by up to 12 km s~^ . 
The lines in the Hermes spectrum also appear deeper, which is 
not simply a consequence of the somewhat lower resolution of 
the UVES spectrum (50 000, compared to 85 000 for the Hermes 
spectrum). Many lines in the UVES spectrum, including the Tc 
lines, seem to suff'er from the often cited "line weakening" effect 
in Mka variables, which causes many spectral lines to become 
much weaker at times, or to disappear completely (Merrill et 



al. 1962 1. This seems to have happened with the near-maximum 
UVES spectrum of R Hya, which we have to re-classify as a 
Tc-rich star This can also be seen as a warning that intrinsically 
weak lines (the Tc lines of R Hya are already at the weak end 
of the Tc line strengths observed in the current sample) may go 
entirely undetected if a Mira star is observed at a pulsation phase 
with severe line weakening. There was no error in the analysis 
of Uttenthale r & Lebzelter| (|2010 ): revisiting the spectmm now, 
we would still classify it as Tc-poor. In view of this situation, we 
scrutinized the present sample to assess whether we could have 
mistaken a Tc-rich star for a Tc-poor one, but conclude that this 
is not the case. 

In Fig. |3] the spectra of three sample stars around the Tc line 
at 423.819 nm are shown. The three stars are chosen to delineate 
a sequence of increasing Tc line strength from absent in T UMi, 
weak in R Hya, and strong in S Her. The lines in R Hya are 
generally weaker than in the other two stars, again illustrating 
the difficulty of detecting the Tc lines in this star. A few of the 
surrounding lines are identified, with identifications being taken 
from |Davis| ( fT947l l. 

The results of our search for Tc are listed in column 5 of 
Table |2] Column 6 summarises what was previously reported 
in the literature about the Tc content of the sample stars. This 
shows that the Tc content of most stars in our sample was un- 
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clear, with there even being contradictory classifications. As dis- 
cussed above, this reflects that the search for the weak Tc lines 
in crowded, low-SNR spectra of red giant stars is not a simple 
and straight-forward task. 

3.1.3. Others-process indicators 

The spectrum of W Dra has basically no signal at the wavelength 
of the Tc lines, hence we searched its spectrum for other indica- 
tors of s-process enrichment. The bands of the ZrO molecule 
are such an indicator, because Zr is also an s-process element. 
The formation of ZrO is also facilitated by an enhanced C/O ra- 
tio, because at C/O ^ 1 the molecular equilibrium is shifted 
to ZrO and LaO, at the expense of TiO and VO. However, in 
contrast to Tc, Zr is a stable element, hence is not an indica- 
tor of a recent s-process and 3DUP. An enhanced Zr abundance 
can also be caused e.g. by the mass transfer of s-process en- 
riched material from a binary companion, which would now be 
a white dwarf. The mass transfer could have happened billions 
of years ago. While Tc would have decayed by now, Zr would 
still be present in the atmosphere of the now evolved secondary 
star. This binary mass transfer scenario is the probable mecha- 
nism leading to Tc-poor S-type stars (Van Eck & Jorissen 1999 1. 
The absence of Zr enrichment (and hence of ZrO bands) would 
however indicate that neither intrinsic nor extrinsic s-process en- 
richment has taken place in a given star. 

Particularly strong and weakly blended bandheads of ZrO 
are to be found in cool star spectra at wavelengths of 647.4 and 
649.5 nm. The combined spectrum of W Dra has enough flux at 
this wavelength to detect the ZrO bands if they were present, but 
apparently they are not (Fig. |4]). In addition, the spectra of all 
other sample stars were inspected for the ZrO bands, but only 
S Her definitely shows therrjj As seen from Fig. |2j S Her is also 
the star with the strongest Tc lines in our sample, so it is clearly 
of MS spectral type. Since ZrO is not detectable in three of the 
Tc- and O-iich sample stars, the absence of ZrO bands does not 
imply the absence of Tc. With ZrO being absent in W Dra, we 
cannot draw a firm conclusion about the Tc content of this star 
We can only conclude that W Dra is an O-rich, M-type star, with- 
out strong s-process (and C/O) enhancement. 



3.1.4. Litliium abundance 

In contrast to Tc, we can determine for Li not only its presence 
but also its abundance, because state-of-the-art model spectra of 
cool giants are able to reproduce observed spectra much more 
closely in the region of the Li line than around the Tc lines. The 
most important stellar parameter that determines the strength of 
the Li I resonance doublet at 670.8 nm is the effective tempera- 
ture. For the 0-iich stars in our sample, the temperature can be 
determined from the strength of the TiO bands. Our method for 
temperature determination involves a grid of COMARCS model 
atmospheres and spectra ( Aringer et al. 2009| l, and closely fol- 
lows the approach of Garcia-Hernandez et al.'('2007). Details of 
the method can be found inlUttenthaler & Lebzelter (2010), and 
here we only outline the main steps of the procedure. 

A grid of synthetic spectra based on model atmospheres with 
temperatures between 2600 and 3600 K was calculated. The step 
width in temperature was lOOK, and two models were added at 
2950 and 3150K. A number of other parameters, which have a 
much smaller impact on the spectrum than the temperature, have 
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As discussed in Sect, 



3.1.1 BH Cru is C-rich, which leaves no trace 



of O-bearing molecules such as ZrO. 



Fig. 4. Spectra of S Her (top) and W Dra (bottom) around the 
ZrO bandheads at 647.4 and 649.5 nm, smoothed by a nine-point 
running box-car The theoretical wavelengths of the bandheads 
are marked by dashed lines. The feature at 647.9 nm is a TiO 
bandhead. 

been fixed beforehand to realistic values. The adopted parameter 
values are: M = IMq, log^[cms-2] = 0.0, Z ^ ^ 0.0156, 
C/O=0.48 (solaiQ), and micro-turbulence ^ - 2.5 km s '. The 
spectra were calculated for regions that contain band heads of the 
TiO molecule, namely 668.4 - 674.0 (which also includes the Li 
line itself) and 701 - 710nm. Before fitting the synthetic spectra 
to the observed ones, they were smoothed with a Gaussian ker- 
nel to the resolution of the respective instrument. An additional 
macro-turbulence had to be adopted, with the maximum value 
of 9 km s ' for R Nor. The synthetic spectra were then fit to the 
observed spectra using ?ix^ minimisation method. The tempera- 
ture of the model spectrum that most closely reproduced a given 
observed spectrum was adopted as the effective temperature of 
that star. These temperatures are listed in column 7 of Table |2] 

Three exceptions to this general procedure have to be noted. 
The Coralie spectrum of R Nor does not reach beyond 690 nm, 
thus the temperature was determined only from the region 668.4 
- 674.0 nm. RU Vul has too weak TiO bands to reliably deter- 
mine the temperature from them. From its mean spectral type 
(M3) and the calibration in |Fluks et al.| ( |1994| ), a temperature of 
3666 K can be inferred. We adopted 3700 K for the spectral syn- 
thesis. At this temperature, the TiO bands predicted by the syn- 
thetic spectra (which are based on solar metallicity) are much 
too strong. However, the temperature may also not be much 
higher than this because the TiO bands, the classification cri- 
terion of M stars, are clearly present in RU Vul, albeit weak. 
The observed spectrum in the two spectral pieces can only be 
fit reasonably well if we assume a metallicity 1 .5 dex below 
the solar value. A comparison with the hotter, metal-poor gi- 
ant Arcturus (Feff = 4290 K, [M/H] = -0.5) also shows that 
RU Vul has much weaker metal lines. Although we did not make 
a full metallicity determination, we suggest that RU Vul is a quite 
metal-poor star. Except for RU Vul, the temperature range of 
our sample stars is fairly narrow. The C-star BH Cru does not 
exhibit the temperature-sensitive TiO bands, so we adopted the 
temperature of 3000 K reported by |Zij Istra et al . | ( j2004) . For the 

- A C/O ratio fixed to the solar value can be problematic for MS- 
or S-stars. However, we have only one MS-star in our sample (S Her, 
which shows ZrO bands), and it turns out that its Li line is undetectable, 
independent of the assumed C/O ratio. 
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Table 2. Results of our analysis. 



Name 


Change 


Sp. 




Tc 


Tc lit. 




loge(Li) 


d\nP/dt 


P 


P 


J-K 


Mk 






type 


type 


(kms-') 






(K) 




(10-^yr-') 


(d, T05) 


(d) 








(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


BH Cru 


SCh 


C 


-3.7 


yes 




3000 


-0.7 


+3.70 


530 


524.4 


1.644 


-8T69 


-5:'59 


RU Vul 


SCh 


M2-4 


-67.9 


no 




3700 


<-1.5 


-7.11 


112 


108.8 


1.223 


-6^88 


-3?89 


TUMi 


SCh 


M4-6 


+8.7 


no 


ul 


3300 


<0.0 


-8.47 


240 


229.1 


1.553 


-7;'44 

-miz; o 


-4?36 



dbfl^ 
no' 

RHya CCh M6-9 -15.9 yes yes* 3100 <0.0 -0.71 385 376.6 1.261 -8?16 -5;H7 

prob" 
no^ 

? ... 2950 <0.0 +1.03 282 289.6 1.340 -7;'78 -4';'74 

^10 ~. 3200 +4:6 +029 5T0 496:2 OTO ^8?55 -5":' 50 

yes yes^ 3300 <0.0 +0.06 307.5 304.1 1.166 -7;'85 -47'89 



WDra 


CCh 


M3-4 


-3.6 


RNor 


MCh 


M3-7 


-28.5 


SHer 


MCh 


M4-7.5S 


-17.5 


SOri 


MCh 


M5-10 


+34.5 


T Cep 


MCh 


M5-9 


-15.7 


WHya 


MCh 


M7-9 


+40.5 


THer 


const. 


M2-8 


-91.9 



yes dbfl2 2950 +0.8 +0.35 430 433.4 1.352 -8?39 -5?49 

yes prob' 3000 +0.5 +0.12 390.9 386.6 1.328 -87'19 -5?16 

yes^ 

no prob' 3300 +0.5 ... 361 388.6 1.459 -8?20 -^25 

no* 

no dbfl^ 3400 +2.0 +0.04 164.9 163.9 1.222 -6^'96 -4';'20 

Note: The meaning of the columns is as follows: (1) Name listed in Simbad; (2) Period change type; (3) Range of spectral type; (4) Heliocentric 
radial velocity; (5) Presence of Tc, this paper; (6) Presence of Tc, previous studies; (7) Effective temperature; (8) Abundance of Li; (9) Rate 
of period change; (10) Period as given in |Templeton et ar] ( |2005^ ; (11) Period found in this study; (12) Mean J-K colour; (13) Absolute K 
magnitude; (14) Bolometric magnitude. References of previous studies on Tc content: l:|Little-Marenin & Littiel ( |1979^ , 2: |Little et al.H1987l l, 3: 
[Vanture et aLl ( [T99T) , 4: |Merrill| ( [T952) , 5: [Uttenthaler & Lebzelter| ( [20T0) , 6: |Lebzelter & Hron| ( [1999| l. Abbreviations mean: "prob" = probably, 
"dbfl" = doubtfiiL 



spectral synthesis, we used model no. 443 from Aringer et al. 
(|2009j, which has T^ff = 3000 K, a C/O ratio of 1.05, and other- 
wise identical parameters to our grid models. A C/O ratio higher 
than 1 .05 does not significantly change the spectral appearance 
in the vicinity of the Li line, hence has very little influence on 
the measured Li abundance. 

In the next step, spectra with varying amounts of Li were 
synthesised based on the best-fit model atmosphere. The hyper- 
fine structure of the Li line was taken into account in the spectral 
synthesis, but not the isotopic shift (i.e. only the contribution by 
^Li was included). This effect can be neglected here, because the 
*Li abundance can be expected to be very small, and other ef- 
fects such as the (velocity) structure of the atmosphere, blending 
with TiO lines, etc., play a much more important role. The Li 
abundance that yielded a minimum in the residuals of observed 
- calculated flux was adopted as the measured Li abundance and 
is reported in column 8 of Table |2] The uncertainty in the Li 



abundance i s however large, up to 0.6 dex ( Garcia-Hernandez 
etaTl [2007} [Uttenthaler & Lebzeltef| '2010'). The largest part of 



this uncertainty is based on the uncertainty in the temperature 
(0.4 dex), the other stellar parameters and flaws in the (TiO) line 
lists contribute to the rest. 

Because of the large opacity of the TiO molecule relative to 
that of C2 and CN at the wavelength of the Li doublet (a factor 
of ~10 difference at Tgff - 3000 K!), the Li detection thresh- 
old for O-rich stars is much higher than for C-rich stars. Hence, 
the abundance of log e(Li) = -0.7 measured in BH Cru would 
not be detectable in M-type stars of similar temperature. Kipper] 
[& W allerstein] ( |1990| l measure a Li abundance in BH Cru of 
-1.5 dex relative to the Sun, or -0.4 on the log e scale, in fairly 
good agreement with our result. The good agreement is remark- 
able because BH Cru has changed its spectral type considerably 
over the past few decades, and our spectrum was taken almost 
26 years after the one used by Kipper & Wallerstein (1990i. 
Furthermore, the synthetic spectrum based on the model cho- 



sen for BH Cru (C/O = 1.05) closely reproduces the molecular 
lines (C2, CN) in the observed spectrum in the vicinity of the Li 
line, which is additional support for a C/O ratio > 1 in that star 

Irrespective of the metallicity of RU Vul, the Li line in this 
star is undetectable. Because of its higher temperature and less 
crowded spectrum, an upper limit of loge(Li) < -1.5 appears 
reasonable. 



The only period change class with noteworthy amounts of 
Li are the meandering Miras. In this group, only S Her does not 
appear to have Li, whereas R Nor has a huge amount of it. The 
measured abundance is log e(Li) = +4.6, which is higher than 
the solar photospheric abundance by a factor of ~3200! The ob- 
served spectrum of R Nor together with a synthetic spectrum to 
fit the Li line is shown in Fig. |5] With the adopted abundance, 
only the red wing of the line is reproduced well, whereas the 
blue wing would require an even higher abundance. However, 
some of the absorption in the blue wing probably stems from 
circumstellar absorption in the expanding stellar wind. This is 
supported by the profile of the Nal D doublet, which also ex- 
hibits a strong blue-shifted absorption component. We did not 
attempt to fit the blue-shifted component in the Li line, because 
our model atmospheres do not include circumstellar shells. The 
very high Li abundance in R Nor suggests that it is undergoing 
HBB, and that it is a faMy massive star. R Nor is in some re- 
spects similar to the very Li-rich S-star V441 Cyg, but is not 
s-process enriched (Uttenthaler & Lebzelter 2010| l. 



Interestingly, the Tc-poor, constant-period Mira T Her also 
contains a considerable amount of Li. It is thus similar to the 



Galactic bulge AGB star Plaut 3-45 discussed in Uttenthaler et 
al. ( |2007b i. This star could be another example of a low-mass 
AGB star experiencing extra-mixing, thereby producing Li. 
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Fig. 5. Observed spectrum of R Nor around the Li line doublet 
(black graph), together with a synthetic spectrum (grey graph). 
The synthetic spectrum is based on a COMARCS model atmo- 
sphere with Teif = 3200 K, log g - 0.0, solar metallicity, and so- 
lar C/O ratio. A Li abundance of log e(Li) - +4.6 was assumed 
in the spectral synthesis. The dotted vertical lines indicate the 
laboratory wavelengths of the Li line doublet. All other features 
in the plotted range are due to TiO. 



3.2. Present-day periods 

Most of our sample stars have experienced a considerable 
change in pulsation period in the past few decades and cen- 
turies. To check whether or not this period change has contin- 
ued in recent years, we analysed the AAVSO visual light curves 
to determine the present-day periods. We downloaded from the 
AAVSO web pag^the visual and V-band observations between 
JD 2 452 000 and JD 2 455 500, i.e. between 3 1 March 2001 and 
3 1 October 2010, and completed a Fourier analysis on the unpro- 
cessed data with the program period04 (Lenz & Breger 2005| l. 
The detected periods (in days) are summarised in column 1 1 of 
Table |2] The peaks in the Fourier spectrum are fairly sharp, the 
Monte Carlo simulation performed with period04 yielded period 
uncertainties of or less than 1 day. Only the Fourier spectrum of 
RU Vul has no clearly highest peak. The one at 108.8 d is only 
barely higher than other peaks, and is broad and asymmetric. 
We thus conclude that no clear periodic pulsation is present any 
more in this star. 

[Templet on et al. | (j2005|l, who use AAVSO data up to about 
1 1 September 2003, do not list present-day periods, but instead 
list average periods over the entire AAVSO time series avail- 
able for a particular star (which can be between a few decades 
and a century long). Nevertheless, they do provide diagrams of 
the period evolution for the stars with the largest period change, 
which are useful to check whether the period change has contin- 
ued. The last period value plotted in these diagrams is the one 
listed in column 10 of Table [2] The meandering Miras S Her and 
T Cep and the "constant" Mira T Her were not found to have 
a large period change (on average over the total time span of 
AAVSO data) by Templeton et al. ( 2005 1, thus no such plots are 



available for them, and we instead list the average period given 
Templeton et al. ( 2005|l. W Hya was not included in the sur- 
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Fig. 6. AAVSO light curve of T UMi since JD 2452000. The 
data points are 10-day means. The previously dominant period 
of =s 230 d has disappeared in the most recent data. 



discuss it further in Sect. 4. 1 In any case, we decided to list in 
Table [2] the period of W Hya as given in the General Catalogue 
of Variable Stars ( Samus et al.[|2009] l. In this table (column 9), 
we also indicate the rate of period change, d In P/dt, defined as 
the slope of the period divided by the average period (Templeton 
et al.) 2005) . For RU Vul, this rate is taken from Templeton et al. 
( 2008[ l, and for W Hya the rate is missing. 

From these data, we can make the following observations. 
Among the SCh change type stars, BH Cru seems to have 
stopped or maybe even reversed its period increase. A decel- 



eration of the period increase was proposed by jTempleton et aL 



emp 



(2005 1. RU Vul has continued its period decrease (also found by 



eton et al. 2008 1, and the pulsation is very weak now. It 
is not detectable any more by visual observations in the most re- 
cent data. The most dramatic evolution in recent years has how- 
ever been seen in T UMi. In the time span investigated here, 
its dominant period decreased to 229 d and did not stop there. 
In the most recent data we see that the star might even have 
switched its pulsation mode from Mira-like to semi-regular vari- 
ability around JD 2 454 600 (May 2008), as the dominant period 
made a jump to 113.6 d! A period ratio close to 2 can be ex- 
pected for a switch from fundamental mode (Mira-like) to first 
overtone (SRV) pulsation (Wood et al. 1991 ). At the same time, 
the full amplitude decreased to about 1™5. According to the clas- 
sical (but arbitrary) definition of AV > 2™5, this small amplitude 
does not qualify T UMi as a Mira any more, but as SRV. The 
light curve of T UMi in the time span investigated here is dis- 
played in Fig.|6] Continued photometric monitoring will be very 
important to follow the exciting evolution of this star. 

Among the CCh Miras, R Hya has continued to decrease its 
period, and W Dra has continued to increase it, while R Aql's 
period remained practically constant. This is somewhat in con- 
tradiction to the luminosity evolution predicted for these stars 
by |Wood & Zarro| ( |1981) , as the period of W Dra is expected 
to stop increasing at some point, whereas R Hya and R Aql are 
expected to continue their decrease, albeit at a somewhat slower 
rate. The development on these short timescales may however 
not be regarded as a stringent constraint of the TP model em- 



ployed by Wood&Zarro(1981i. Among the MCh Miras, R Nor 



and W Hya are the stars that display the strongest changes. The 
former decreased its period in the recent years, while the latter 
increased it considerably. 



in 



vey of Templeton et al.'(2005), probably because it is classified 3 3 Luminosities 



as SRV in Simbad. This classification is controversial, and we 



http://www.aavso.org 



To compare the evolutionary states in the HRD of the dififerent 
period change classes, as well as of the Tc-rich and the Tc-poor 
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stars, we determined the absolute bolometric magnitudes (lumi- 
nosities) of our sample stars. In principle, this could be done 
very easily using a period-luminosity relation for Mira stars, 
e.g. the one from Guandal ini & Busso| ([2008 their Fig. A.l). 
However, our sample stars do not necessarily closely follow this 
relation, because of the change in the pulsation period. For ex- 
ample, |ZyIsttaetir('2007) argue that the near-IR data of BH Cru 
is consistent with evolution at constant luminosity, and that any 
increase in luminosity did not exceed 10 percent. If this were 
true, then the star would now be fainter than other stars of a sim- 
ilar pulsation period that do follow the Mira period-luminosity 
relation. The only way to determine luminosities independently 
of the pulsation period is by the measurement of the paral- 
lax of the stars. Unfortunately, parallaxes of high significance 
(cTji/tt < 0.5) are available for only four of our sample stars 
dPourbaix et aLi[2003l ). 

In view of this situation, we decided anyway to infer the lu- 
minosity from the pulsation period. In a first step, we used the 
log P — K relations of |Riebel et al. ( |2010[ ) to determine the abso- 
lute K magnitudes. As these relations are based on LMC stars, 
we assume that the same relation is applicable to Galactic AGB 
stars. This will likely introduce only a small uncertainty in the 
final luminosity. For the carbon star BH Cru, the relation for se- 
quence 1 (classical Mira sequence) of C-rich stars, and for all 
other stars sequence 1 of the O-rich stars was used. As the pe- 
riod of RU Vul is no longer clearly detected in the most recent 
photometric data, we adopted the period of 155 d that it had be- 
fore the period decline ( [Templeton et al.[ |2008| and forced it to 
follow also sequence 1 of the fundamental mode pulsators. This 
choice is justified in Sect. 



4.1 



A distance modulus to the LMC 
of 18™5 was assumed. To get from absolute /T-magnitudes to the 
absolute bolometric magnitude Mboi, we applied the BC(K) re- 
lations of |Kerschbaum etaL] ( |20101 l based on the J — K colour. 
Quasi-simultaneous photometry in the / and K bands was col- 
lected f rom the works of [Ca tchpole et al. ( 1979[), Fouque et al. 
( fT992l l, [Kerschbaum & Hron. (1994|),"|Whitelock et al.| ( |2006| l 



and the 2MASS catalogue ( Skrutskie et al. 2006| l, transformed 
to the 2MASS photometric system with the relations given by 
Carpenter ( 200 l| l, and averaged to a final value. For some stars, 
only the 2MASS data were available, but for most stars this re- 
sulted in a cycle-averaged J - K colour As in Kerschbaum et al. 
( 2010| l, no correction for interstellar reddening was applied, be- 
cause E(J- K) is only a few 0'"01 for these solar neighbourhood 
stars and variability and individual circumstellar extinction have 
a much larger influence. The BC(K) correction applied ranged 
between 2™76 and 3™ 11. The mean J - K colour, the absolute K 
magnitude, and the resulting bolometric magnitude Mboi are re- 
ported in columns 12, 13, and 14 of Table [2] The Mboi derived in 
this way differ on average by less than 0"08 from what is found 
by using the relation of jGuandalini & Busso| ( |200"8] ). We estimate 
a typical uncertainty of ~(J^16 in the bolometric magnitude from 
combining the uncertainties in the log P-K and BC(K) relations. 
The variable period of the stars might introduce an additional un- 
certainty, which is hard to quantify. 

For the four stars in our sample with a good parallax (R Aql, 
R Hya, T Cep, and W Hya), we can compare the absolute K- 
magnitude derived from the log P-K relation with that derived 
from the parallax. Only in the case of W Hya do the values not 
agree within the error bars, from the log P-K relation being 
brighter by (y?36 than the brightest Mk allowed by the Icr uncer- 
tainty in the parallax. The average difference for the other three 
stars is 0.^18, in good agreement with the error estimated for the 
luminosities adopted here. For R Nor, the uncertainty in the par- 
allax is slightly larger than half the parallax itself. Nevertheless, 




Fig. 7. HRD with our sample stars. The symbols are as follows: 
Triangles represent the SCh stars, squares the CCh Miras, circles 
the MCh Miras, and the rectangle represents the constant Mira 
T Her. Open symbols are the Tc-poor stars, filled symbols Tc- 
rich. The gray square is W Dra, for which no definite conclusion 
on its Tc content could be drawn. Stars discussed individually in 
Sect.l4.1lare identified with their names. 



a comparison of the Mk magnitudes derived with the two meth- 
ods is interesting and suggests that the true parallax is probably 
near the lower limit of the Icr range, hence the distance is prob- 
ably near the upper limit of 800 pc (see also Sect. |4.1| i. 

The resulting HRD of our sample stars, using J - K as a 
proxy for the temperature, is shown in Fig. |7] As already noted 
by Uttenthaler et al.| ( 2007a i, the Tc-rich stars seem to form the 
upper envelope of the distribution. They are the brightest stars 
at any given J - K colour (or temperature). The only exception 
here is R Nor, a very bright but definitely Tc-poor star This is 
further independent evidence that R Nor is a massive star with 
HBB (see Sect.|4T|i. From its location in the HRD, W Dra falls 
among the Tc-poor stars, which suggests that it is also Tc-poor. 



4. Discussion 

Before we discuss the evolutionary states of the different period 
change classes, we first comment on the evolutionary states of a 
few remarkable individual objects in the sample. 

4.1. The evolutionary states of individual objects 

BH Cru displays Tc in its spectrum, a clear sign that it is on 
the TP- AGB. An s-process enhancement in that star has already 
been reported by'Abia & Wallerstein ( 1998i The pulsation pe- 
riod of BH Cru has increased considerably since 1970, but has 
been stable for some 15 years, or even slightly decreased. We 



presented in Sections 3.1.1 and 



3.1.4 evidence that BH Cru is a 



carbon star. It was initially classified as being of SC spectral type 



(with C/O essentially unity; Catchpole & Feast 1971[ Keenan 



I& Boeshaar [1980), but its type changed to CS and it experi- 
enced a marked period increase (Lloyd Evans 1985|. [Zijlstra et| 
al.| ( |2004| interpret this spectral change as being due to a de- 
crease in the effective temperature by some 200 K, rather than 
a change in C/O ratio (i.e. by a recent dredge-up of C). The 
same authors obtained a low-resolution spectrum of BH Cru in 
the range 600 - 700 nm and noted that C2 bands are not clearly 
present. However, they also stated that |Loidl et al.| ( [2()01[ ) found 
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that C2 bands are present longward of 700 nm. We compared 
the spectrum of BH Cru to Hermes spectra of well-known car- 
bon stars and found that RS Cyg, a carbon star of spectral type 
C8.2, has an almost identical spectrum. The measured C/O ra- 



tios of C-stars are clearly in excess of 1 (e.g. Lambert et al 
1986 Ledereretal. 2009 [ l. We favour the interpretation that the 



change in spectral type observed in the past few decades is re 
ally due to a recent 3DUP event, and not simply a temperature 
decrease. If that were the case, then every normal carbon star 
could become either a CS or an SC star simply by increasing its 
effective temperature. Dedicated modelling of the atmospheres 
and spectra of cool giant (Mira) stars with C/O ratio close to 
unity, including non-equilibrium chemistry, is required to gain 
more certainty in this respect. We note that an increase in the 
C/O ratio from below to above unity at constant luminosity nat- 
urally leads to an increase in the pulsation period, because of 
the associated change in chemistry and opacities. This is sug- 
gested by the linear pulsation models employed by Lebzelter & 
|Wood| ( |2007| l, who take into account the effect of the C/O ratio 
on the opacities and the atmospheric structure of pulsating AGB 
stars. This effect remains to be confirmed with more realistic 
non-linear pulsation models. If confirmed, this would mean that 
the observed period change in BH Cru can actually be caused 
by the dredge-up of carbon and the thereby changing opacities, 
rather than by the associated TP itself. 

RU Vul no longer displays any clearly detectable pulsation, 
and shows signs of neither Tc nor Li. We find that it has a very 
low metalhcity of [M/H] ^ -1.5. This agrees with the finding 
of Mennessier et al. (2001 1, who assign it to the extended (thick) 



disk, which is an old, metal-poor population. The same authors 
determine an absolute IRAS 12 magnitude of RU Vul (based on 
the statistical properties of the extended disk sub-sample), from 
which we infer a distance of 2070 pc to that star. If we assume a 
fundamental mode pulsation with a period of 155 d (Templetqn] 
let al.||2008] l, use relation 1 of Rie bel etal7|p010 ) to derive an ab- 
solute /T-magnitude, and combine this value with the observed 
2MASS ^T-magnitude, we arrive at exactly the same distance. 
This provides good confidence that assigning RU Vul to the fun- 
damental mode pulsators is the correct choice. |Kerschbaum & 



[Hron ( 1992) find that the SRa's are not a distinct class of vari 
ables but a mixture of 'intrinsic' Miras and SRb's. Before its 
period and amplitude decline, RU Vul exhibited a quite regular 
light change with a visual amplitude of about 2™0, somewhat be- 
low the (arbitrary) definition limit of Miras, and was classified 
as SRa. Furthermore, many short-period, small-amplitude Miras 



are known to belong to an old, metal-poor population (Hron 
1991[ l. We conclude that RU Vul is also representative of this 



very same population, but it failed to be classified as Mira be- 
cause of its too low pulsation amplitude. Its supposedly high age 
also suggests that RU Vul is of relatively low mass (~ IMq). 

T UMi is an oxygen-rich, M-type star without any signs of 
s-process enrichment. It has thus not undergone a 3DUP event. 
The change in the pulsation behaviour of this star is dramatic. 
Before 1980, the period was fairly stable at ~315d, when it 
started to decline at a high rate. By mid 2008, its period had de- 
creased to ~230d, when it seems to have switched the pulsation 
mode from Mira-like to semi-regular variability (Fig. |6|l. The 
present-day period is around 113 d, with a visual amplitude be- 

1999[ l already spec- 



low the definition limit of Miras. iWhitelock 



ulate that T UMi could drop out of the Mira instability strip in the 
HRD, which we confirm herewith. [Mattel & Foster] ( |l 99 5 ) sug- 
gest that it is in the early helium-shell flash phase where the pe- 
riod is expected to change most rapidly as the luminosity drops. 
An indication for a luminosity drop in T UMi was found in the 



detailed light curve analysis by Szatmary et al. ( 2003| l, who also 
suggest that the period decrease should stop in the near future, 
which is however not (yet) the case. However, from the lack of 
Tc (and also Li) it seems that T UMi is not as highly evolved as 
other stars in the present sample, or that it is not massive enough 
to undergo a 3DUP event. If a recent TP is the cause of its period 
change, it is possibly an early (weak) TP. 

R Hya is the longest-known Mira to have a steadily decreas- 
ing pulsation period. At the time of its discovery around AD 
1700, it had a period of ~500d, which then decreased to its 
present value of ~376 d, and might well continue to decline in 
the future. We find Tc at a low level in this star, which confirms 
its evolutionary state on the TP- AGB. We thus conclude that the 
TP scenario for the period change in this star is likely. 

The high Li abundance in R Nor (Sect. [3.1. 4| l could be a 
result of hot bottom burning (HBB) going on in this star This 
would require the star to have a mass of ~ 4Mq. (Lebzelter et| 
[aL] ( |2005j ) suggest that R Nor is relatively massive (3 - 5Mo), 
because it has a secondary maximum in its light curve, which is 



also found in luminous LMC Miras (see also McSaveney et al. 



2007). This would be in line with our finding. In the LMC, the 
very Li-rich AGB stars are mainly found in the luminosity range 
-6^0 > Mboi > -7^2 ( [Smith et al.| [T995| ). R Nor is somewhat 
below this limit (Mboi = -5"50, Sect. |4.1[ ), with the uncertainty 
that we have to rely on the period to determine its luminosity. 
Another indicator of a high mass can be a relatively small dis- 
tance from the Galactic plane. The Galactic latitude of R Nor is 
5° and its parallax n = 2.76 ±1.71 mas (Pourbaix et al.||2003|l. 



Even though the parallax is very uncertain, this would place 
R Nor at a distance from the Galactic plane between 20 and 70 pc 
(where the larger distance is the more likely one, see Setc. 3.3 1. 



For Miras with periods in the range 300 - 400 d, the exponential 
scale height is close to 240 pc Pura & Kleinmann 1992| l, hence 
R Nor is relatively close to the Galactic plane. Furthermore, 
Garcfa-Hernandez et al. (2007) find no s-process enrichment in 
their sample of Li-rich, massive Galactic AGB stars. This could 
be explained by the presence of a massive convective envelope, 
which would strongly dilute any material dredged-up from re- 
gions close to the core. In addition, 3DUP could be very inef- 
ficient in high mass stars. An additional piece of evidence of 
a high mass would be an enhanced rubidium (Rb) abundance, 
which is synthesised by the ^^Ne(a,n)^^Mg neutron source reac- 



tion (Garcfa-Hernandez et al. 2006 1. Unfortunately, our Coralie 



spectrum of R Nor does not cover the Rb line at 780 nm, so we 
cannot measure its abundance in this star. Nevertheless, there is 
compelling evidence that R Nor is of high mass, probably around 
4 Mo. 

Whether W Hya is an SRV or a Mira is not entirely clear. 
This question was discussed in detail by Nowotny et al. (20T0|. 
From its light variation and location in the P-L diagram, W Hya 
would indeed be classified as a Mira. From the velocity varia- 
tions in its atmosphere, on the other hand, it can clearly be as- 

(!2005) conclude 



signed to the group of SRVs. Lebzelter et al. 



that a high mass might be the reason for the remarkable veloc- 
ity behaviour. Although some Li is present in the atmosphere of 
W Hya, the mass of this star is probably too low for efficient 
operation of HBB. 

4.2. The evolutionary states of the period change classes 

The SCh group in our sample is very heterogeneous: The period 
of the C- and Tc-rich Mira BH Cru is (or has been) increasing, 
the O-rich, Tc-poor stars T UMi and RU Vul have decreased 
their periods considerably. In addition, the derived luminosities 
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(Sect. [33] l differ widely. These three stars differ quite fundamen- 
tally from each other, including in terms of their evolutionary 
states. Nevertheless, they show similarly large period change 
rates (in absolute terms). If their period change is caused by a 
TP, the inhomogeneity of the presence of Tc in the SCh group 
suggests that a TP with an associated large period change can 
happen before a star undergoes a 3DUP event. 

From the viewpoint of chemistry and luminosity, the CCh 
and the MCh groups are much more homogeneous than the SCh 
group. All members of these two groups studied here are late 
M-type stars. However, both the CCh and the MCh groups are 
inhomogeneous in terms of the presence of Tc, even though for 
W Dra no definite conclusion could be drawn. There is some ev- 
idence that the MCh group of stars is on average more massive 
than the other two groups. The periods in this group are fairly 
long, between 300 and 500 d, although this is in some contra- 
diction to |Zijlstra & Bedding ( 2002 1, who note that all mean- 
dering Miras have periods of 400 d or longer Four of five stars 
in this group have at least some Li, and three display Tc. The 
presence of Tc certifies that they are definitely in the TP-AGB 
phase. We do, however, not suggest that the period change in the 
MCh group is connected to a recent TP, as their periods change 
back and forth on a much too short timescale. One of the alter- 
native explanations presented in the introduction (feedback be- 
tween pulsation and atmospheric structure or chaotic behaviour) 
might hold for them. 

5. Summary and conclusions 

From our search for Tc and other s-process indicators, as well 
as for Li, in the spectra of Mira variables with changing pul- 
sation periods, we can draw several conclusions. In our sample 
of twelve stars, have identified five that display Tc absorption 
lines, indicating that they have undergone a 3DUP event follow- 
ing a TP. The most important conclusion, however, is that there 
is no clear correlation between the presence ofTc and the type of 
period change, even in this limited sample. We suggest that the 
two Tc-rich, non-meandering Miras BH Cru and R Hya are the 
most likely to have experiences a recent TP, although their pe- 
riod change rates are quite different. Furthermore, we presented 
evidence that BH Cru changed its spectral type from SC to C, 
which indicates that there has been a recent 3DUP that raised its 
C/0 ratio above unity. There is also accumulating evidence (high 
Li abundance, high luminosity, pulsation behaviour, proximity to 
the Galactic plane) that R Nor is a fairly massive star of ~ 4Mq. 
On the basis of an analysis of AAVSO data of T UMi, we report 
that this star probably switched its pulsation mode from Mira- 
like to SRV very recently. Finally, we suggest that RU Vul is a 
representative of the old, metal-poor thick disk of a consequently 
low mass. 

Continuous photometric monitoring of Mira variables, in 
particularly those with large period changes, is mandatory to 
gain further insight into the AGB phase of stellar evolution. The 
archives of amateur observations such as those at AAVSO are 
an invaluable resource of data for this study. In parallel to pho- 
tometric monitoring, we also recommend spectroscopic moni- 
toring of the stars with changing pulsation period, albeit with 
sparser sampling of a few years. This will be of high importance 
to check whether the period change is also accompanied by a 
spectral change. 

In addition to the prominent and well-studied objects pre- 
sented here, there are more stars with pronounced period 
changes that deserve more detailed scrutiny in the future. 
LX Cyg ( [Zijlstra et al.| |2004| l is a particularly interesting ex- 



ample, which is classified as SC spectral type, just as BH Cru 
was before its spectrum started to change. Assuming that the 
s-process enrichment in LX Cyg is intrinsic, this star would sim- 
ilarly be placed on the TP-AGB. A period change caused by a 



recent TP seems plausible. Furthermore, TT Cen (Whitelock 



| 1999[ ), S Sex ( [Merchan Benitez & Jurado Vargas[|2000|l, as wel 
as the remarkable objects mentioned in'Temp leton et al.| ( |2005 



R Cen, Z Tau, Y Per) deserve a closer look with high-resolution 
spectroscopy, which we have to postpone to a future project. 
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